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The chiral low-energy constants c_d and ce are constrained by means of accurate ab initio calcu- 
lations of the A — 3 binding energies and, for the first time, of the triton f3 decay. We demonstrate 
that these low-energy observables allow a robust determination of the two undetermined constants, a 
result of the surprising fact that the determination of c_d depends weakly on the short range correla- 
tions in the wave functions. These two- plus three-nucleon interactions, originating in chiral effective 
field theory and constrained by properties of the A = 2 system and the present determination of c_d 
and Ce, are successful in predicting properties of the A = 3, and 4 systems. 

PACS numbers: 21.30.-x,21.45.Ff,23.40.-s 



The fundamental connection between nuclear forces 
and the underlying theory of quantum chromodynamics 
(QCD) remains one of the greatest contemporary the- 
oretical challenges, due to the non-perturbative charac- 
ter of QCD in the low-energy regime relevant to nuclear 
phenomena. However, the last two decades of theoreti- 
cal developments provide us with a bridge to overcome 
this obstacle, in the form of chiral perturbation theory 
(xPT) [T]. The xPT Lagrangian, constructed by inte- 






grating out degrees of freedom of the order of ~ 1 GeV 
and higher (nucleons and pions are thus the only explicit 
degrees of freedom) , is an effective Lagrangian of QCD at 
low energies. As such, it retains all conjectured symme- 
try principles, particularly the approximate chiral sym- 
metry, of the underlying theory. Furthermore, it can be 
organized in terms of a perturbative expansion in positive 
powers of Q/A^ where Q is the generic momentum in the 
nuclear process or the pion mass [IJ. Though the subject 
of an ongoing debate about its validity [2, Jj, the naive 
extension of this expansion to non-perturbative phenom- 
ena provides a practical interface with existing many- 
body techniques, and clearly holds a significant value for 
the study of the properties of QCD at low energy and its 
chiral symmetry. 

The chiral symmetry dictates the operator structure of 
each term of the effective Lagrangian, whereas the cou- 
pling constants (not fixed by the symmetry) carry all the 
information on the integrated-out degrees of freedom. A 
theoretical evaluation of these coefficients, or low-energy 
constants (LECs), is equivalent to solving QCD at low- 
energy. Recent lattice QCD calculations have allowed a 
theoretical estimate of LECs of single- and two-nucleon 
diagrams T, while LECs of diagrams involving more 
than two nucleons are out of the reach of current compu- 
tational resources. Alternatively, the undetermined con- 
stants can be constrained by low-energy experiments. 

The strength of xPT is that the chiral expansion is 



FIG. 1: Contact and one-pion exchange plus contact interac- 
tion (a), and contact MFC (fo) terms of xPT. 



used to derive both nuclear potentials and currents from 
the same Lagrangian. Therefore, the electroweak cur- 
rents in nuclei (which determine reaction rates in pro- 
cesses involving external probes) and the strong interac- 
tion dynamics {ttN scattering, the NN interaction, the 
NNN interaction, etc.) are all based on the same the- 
oretical grounds and rooted in the low-energy limits of 
QCD. In particular, xPT predicts, along with the NN 
interaction at the leading order (LO), a three-nucleon 
(NNN) interaction at the next-to-next-to- leading order 
or N^LO [5J IS], and even a four-nucleon force at the 
fourth order (N^LO) [7]. At the same time, the LO 
nuclear current consists of (the standard) single-nucleon 
terms, while two-body currents, also known as meson- 
exchange currents (MEC) , make their first appearance at 
N^LO [m. Up to N^LO both the NNN potential and 
the current are fully constrained by the parameters defin- 
ing the NN interaction, with the exception of two "new" 
LECs, CD and ce- The latter, ce, appears only in the 
potential as the strength of the NNN contact term [see 
Fig. [l](a)]. On the other hand, cd manifests itself both 
in the contact term part of the NN-ir-N three-nucleon 
interaction of Fig. [T|(a) and in the two-nucleon contact 
vertex with an external probe of the exchange currents 
[see Fig.[l](6)]. 
First attempts to determine cjj and ce have used 



the triton binding energy (b.e.) alongside an addi- 
tional strong observable, either the nd doublet scattering 
length 6 , or the *He b.e. [8]. However, this led to a 
substantial uncertainty in the values of the LECs, due to 
correlations between the '^H b.e. and these observables, 
known respectively as the Phillips and Tjon lines. The 
fine-tuning of these observables is very sensitive to the 
structure of the adopted NNN force. Hence small varia- 
tions of the cutoff, different regularization schemes, miss- 
ing terms of the interaction, etc., tend to produce large 
swings in the extracted of cd and ce- A different ap- 
proach was adopted in Ref. [3]. There, a preferred choice 
for the two LEC's was obtained by complementing the 
constraint on the A = 3 b.e. with a sensitivity study on 
the radius of ''He and on various properties of p-shell nu- 
clei. The same interaction was then successfully used to 
predict the ''He total photo-absorption cross section [TU] . 

The need for a complemental determination of these 
LECs is two-fold. First, it would be desirable to per- 
form such a determination within the A < 3 systems, 
so to suppress any additional many body contribution. 
Second, despite recent progress [TT], the NNN potential 
has been fully worked out only up to N^LO, leaving in- 
consistency in the calculation of the wave functions when 
combined with the NN force at N'^LO. Clearly, it would 
be preferable to adopt an observable with minimal depen- 
dence on the short-range part of the wave function. In 
this respect, the relation (mandated by the chiral symme- 
try of QCD) between electroweak processes and NNN- 
force effects offers venues to achieve these goals. This 
relation was established in the context of effective field 
theory [12l [13j [M] , and manifests itself in xPT via the ap- 
pearance of cd in both the NN-n-N diagram of Fig. [l](a) 
and the one in Fig. [l](6). In particular, Gardestig and 
PhiUips [13], suggested the triton beta-decay as one of 
the electroweak processes that could be used as input to 
fix the strength of the NNN force. It is the purpose of 
this Letter to undertake this task and show that by us- 
ing the triton half life, as well as the A = 3 b.e., one can 
constrain the two undetermined LECs within the three- 
nucleon sector by means of fully converged ab initio cal- 
culations. We demonstrate that this determination is ro- 
bust. The resulting chiral Lagrangian predicts, without 
any free parameters, various ^ = 3, and 4 properties. 

The triton is an unstable nucleus, which undergoes 
/3-decay with a "comparative" half-life of (/T'i/2)t = 
(1129.6 ± 3) s [15]. This quantity can be used to extract 
an empirical value for (E^) = K^ReWE^ 11^11)1 ^ [TT], 
the reduced matrix element of the J = 1 electric multiple 
of the axial vector current, through 
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Here, ii' = 27r^ In 2/TOg (with rUe the electron mass), Gv 
is the weak interaction vector coupling constant (such 




FIG. 2: (Color online.) cd-ce trajectories fitted to reproduce 
and ^He experimental b.e. The dotted lines show the 
region for which \1 - {E^)th/ {E^)emp\ is within the ±0.54% 
error-bars. 



that K/G^ = 6146.6 ± 0.6s T5]), fA/fv = 1.00529 [T7] 
accounts for the small difference in the statistical rate 
function between vector and axial- vector transitions, and 
6c = 0. 13% jl7| is a small correction to the reduced matrix 
element of the Fermi operator, calculated between the 
^ = 3 wave functions (which is 1 for this specific case) due 
to isospin-breaking in the nuclear interaction. One can 
use these values to extract (£'j'^)|emp = 0.6848 ± 0.0011. 

The weak axial current adopted in this work is the 
Nother current built from the axial symmetry of the chi- 
ral Lagrangian up to order N'^LO [19 . At LO this current 
consists of the standard single- nuclcon part, which at low 
momentum transfer is proportional to the Gamow- Teller 
(GT) operator, -E^^Ilo = « .9A(67r)"i/2 cr.,T+ , ^j^ej-g 
ai , are spin and isospin-raising operators of the ith nu- 
cleons, and = 1.2695±0.0029 is the axial constant |20j . 
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For this reason, the quantity 
referred to as "experimental", or "empirical", GT. 

Corrections to the single-nucleon current appear at 
N^LO in the form of MEG and relativistic terms. The 
MEG are formed by a one-(charged)-pion exchange, and 
a contact term. While the relativistic corrections are neg- 
ligible for the triton half life, the MEG have a substan- 
tial influence on this /3-dccay rate. This is a reflection of 
the fact that E^ is a chirally unprotected operator [21] . 
Moreover, the strength of the MEG contact term, usually 
denoted by dj^, is related to cd through: 
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Here, is the nucleon mass, and C3 and C4 are LECs 
of the dimension-two ttN Lagrangian, already part of 
the chiral NN potential at NLO. Therefore, one can use 
{E^)\emp as second constraint for the determination of 
C£) and Ce- 

Following the cd-ce trajectory which reproduces, on 
average, the ^ = 3 b.e., as discussed in Ref [3], here, we 
(i) calculate the ^H and '^He g.s. wave functions by solv- 
ing the Schrodinger equation for three nucleons interact- 
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FIG. 3: (Color online.) The ratio {E^)th/{E^)emp using the 
N^LO iViV potential [22] with and without the local N^LO 
NNN interaction [23], and the axial current with and with- 
out MEG {cd,ce are varied along the averaged trajectory of 
Fig. [2]|. The shaded area is twice the experimental uncer- 
tainty. Also shown: results (without NNN force) for the 
phenomenological AV18 potential (with A — 500 MeV im- 
posed in the current), and for the N'^LO A'^A'' potential of 
Epelbaum et al. (with A = 450, 600 MeV, and a 700 MeV 
spectral- function cutoff in the two-pion exchange). 



ing via the xPT NN potential at N^LO of Ref. [22] and 
the NNN interaction at N^LO [5] in the local form of 
Refs. [H [ini US [25 ; (ii) determine for which cd values 
along the trajectory the calculated reduced matrix ele- 
ment of the Ef operator (at N^LO) reproduces {Ei)emp- 

The present calculations are performed in the frame- 
work of the no-core shell model (NCSM) approach [5] [HJ 
[53] [5S]. This method looks for the eigenvectors of the 
Hamiltonian in the form of expansions over a complete 
set of harmonic oscillator (HO) basis states up to a max- 
imum excitation of N^i^^Ml above the minimum energy 
configuration, where is the HO parameter. Thanks 
to the large model-space size adopted (A^max — 40), 
A — 3 b.e. and reduced matrix element of are con- 
verged to less than 0.05%. Note that the same regulator 
FA{q^) = exp(-g''/A'') is used for both NNN terms of 
the interaction and MEG, a process resulting in a local 
chiral NNN force (for relevant parameters and defini- 
tions see Ref. g3]). The A = 3, 4 calculations of Ref. [23] 
were later confirmed by the results of Ref. [26j . provid- 
ing a benchmark for the local chiral NNN force. The 
MEG utilized in this work were validated against those 
of Park et al. 19J. Finally, we tested the implementation 
of the MEG within the NGSM approach by reproducing 
(within 0.1%) the AV18 [27] results for {E\) obtained 
using the effective-interaction hyper-spherical harmonics 
approach [19 . 

The theory to empirical value ratio for the E^ reduced 
matrix element along the averaged constraint of Fig. [2] 
(which reproduces the A = 3 b.e. to about lOkeV p^) 
is presented in Fig. [3] The 1.08% tolerance band high- 
lighted by the shaded area (obtained by doubling the er- 
ror bar) is mainly due to the uncertainties on {Ef-)\emp 



and gA- Besides the full calculation, which appears as 
a solid line, we report also the results of several tests, 
aimed to analyze the sensitivity of the triton half life to 
NNN force and/or MEG. 

First we note the fundamental importance of the ax- 
ial two-body currents in reaching agreement with experi- 
ment. By suppressing the MEG, in the whole investigated 
cd-ce range, the calculations under-predict {Ef)\emp by 
about 2%. The same, almost constant, behavior is found 
when adding to the single-nucleon current only the long- 
range one-pion-exchange term of the MEG, which cor- 
responds to setting = 0. In this case, the theoreti- 
cal results over-predict {Ef)\emp by ~11%. Only when 
adding the contact part of the MEG, which is related to 
the short range weak correlations of axial character, can 
the half-life reach its experimental value. In particular, 
we find that agreement within ±0.54% of the empirical 
value is obtained for — 0.3 < c^i < — 0.1. The correspond- 
ing ce values lie in the range [—0.220,-0.189]. These 
results are summarized by the dotted lines in Fig. |2] 

In a similar spirit, we now study the effect of the sup- 
pression of the NNN force. If we try to calibrate cd 
to reproduce the measured half-life, we obtain a curve 
in close agreement with the results of the full calcula- 
tion [13] (for completeness we show also the curve cor- 
responding to the suppression of both MEG and NNN 
force). Moreover, a quantitatively similar cu dependence 
of (E^) can be obtained using A — 3 wave functions 
produced by the phenomenological AY 18 NN potential 
(without NNN force). It is therefore clear that the half 
life of triton presents a very weak sensitivity to the NNN 
force, and hence to the strength of the spin-orbit interac- 
tion. Thanks to this feature, which might be unique to 
s-shell nuclei, we are confident that the determination of 
ce) and ce obtained in this way is robust. Incidentally, 
the weak dependence of the half life of triton upon the 
NNN force can also explain the success of recent calcu- 
lations done in a hybrid approach, coined EFT* [T^ . 

As the values of the C3 and C4 LEGs are somewhat 
uncertain (see, e.g., Ref. 0), it is important to assess 
to which extent they would influence the determination 
of Cd from the triton half life. While very sensitive 
to the smallest change in C3, the N^LO fit of the NN 
data of Ref. [22J does not deteriorate dramatically for 
3.4GeV"^< C4 < 5.4GeV-^ [2H]- Fig. [s] shows calcula- 
tions (without NNN force) carried out by setting C4 to 
3.4GeV~^ {ttN value [29]) in the axial current, while the 
A — 3 wave functions are still obtained from the N'^LO 
NN potential of Ref. [H] (where, in GeV"\ C3 = -3.2 
and C4 = 5.4). We find that the use of the lower C4 value 
produces a shift (^0.3) towards more positive cd values. 

Finally, Fig. |3] shows two additional curves (without 
NNN force) obtained using the N^LO NN potential of 
Ref. [3D] with A = 450 and 600 MeV (700 MeV two- 
pion exchange spectral-function cutoff) and the parame- 
ters defining it (particularly, in GeV""'^, C3 = —3.4 and 
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TABLE I: Calculated ^'H, ^He and *He g.s. energies (in MeV) and point-proton radii (in fm), obtained using the N^LO NN 
potential [21] with and without the local N'^LO NNN interaction [23] with cr>~~0.2 and ce = —0.205, compared to experiment. 









^He 






*He 






-Bg.s. 










-7.852(4) 


1.651(5) 


-7.124(4) 


1.847(5) 


-25.39(1) 


1.515(2) 


NN+NNN 


-8.473(4) 


1.605(5) 


-7.727(4) 


1.786(5) 


-28.50(2) 


1.461(2) 


Expt. 


-8.482 


1.60 


-7.718 


1.77 


-28.296 


1.467(13) [H] 



C4 = 3.4). As one could expect, the extracted value of 
the NNN LEG cd depends on the choice of the cutoff. 
However, we observe that the A — 450/700 MeV poten- 
tial gives comparable results as the N'^LO NN potential 
of Ref. [22] and AV18, indicating that the determination 
of C£) depends mainly on the MEG cutoff and weakly on 
the cutoff imposed in the nuclear potential. 

With this calibration of cd and ce, for this potential, 
in principle, any other calculation is a prediction of xPT. 
In Table |T] we present a collection of A = 3 and 4 data, 
obtained with and without inclusion of the NNN force 
for CD — —0.2 [ce = —0.205), a choice in the middle of 
the constrained interval. Besides triton and '^He g.s. en- 
ergies, which are by construction within few keV from 
experiment, the NN + NNN results for the ^He are in 
good agreement with measurement. Note that a particle 
g.s. energy and point-proton radii change minimally with 
respect to variations of cd in the interval [—0.3,-0.1], 
and the results at the extremes are both within the nu- 
merical uncertainties quoted in Table |T] This result is not 
inconsistent with the study of Ref. [9 , which showed pref- 
erence for Cd ^ — 1, since for p-shell nuclei one expects 
the (neglected) higher-order NNN force terms to affect, 
probably through a shift, the value of cd [11] . 

Summarizing, we have used the A = 3 b.e. and the half- 
life of triton to constrain the undetermined N'^LO xPT 
parameters of the NNN force. We have demonstrated 
the robustness of the constraint on ce by showing the 
weak sensitivity of the (E^) matrix element with respect 
to the NNN force. In particular, we find —0.3 < cd < 
-0.1, and, correspondingly, -0.220 < < -0.189. The 
latter is expected to change due to N'^LO terms of the 
NNN interaction, which were not included thus far. In 
conclusion, we have identified a clear path towards de- 
termining the NNN force that, once the NN interaction 
will be pinned down, will pave the way to parameter-free 
predictions of QGD in the consistent approach provided 
by xPT. 
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